DFTUZ /98/10 



EXACT SOLUTION OF A ELECTRON SYSTEM 
COMBINING TWO DIFFERENT t-J MODELS 

J. Abad and M. Rfos^^) 
Departamento de Fisica Teorica, Facultad de Ciencias, 
Universidad de Zaragoza, 50009 Zaragoza, Spain 



Abstract 

A new strongly correlated electron model is presented. This is formed by two types of 
sites: one where double occupancy is forbidden, as in the t-J model, and the other where 
double occupancy is allowed but vacancy is not allowed, as an inverse t-J model. The 
Hamiltonian shows nearest and next-to-nearest neighbour interactions and it is solved by 
means of a modified algebraic nested Bethe Ansatz. The number of sites where vacancy is 
not allowed, may be treated as a new parameter if the model is looked at as a t-J model 
with impurities. The ground and excited states are described in the thermodynamic limit. 
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1. Introduction 



Strongly correlated systems are very interesting in view of their relation with high 
Tc superconductivity. It is very important to study the one-dimensional systems because 
they may share properties with two-dimensional ones The t-J model was proposed 
by F.C. Zhang and T.M. Rice and it describes electrons on a lattice excluding the 
double occupancy of any site, in opposition to the Hubbard model 0, where the double 
occupancy is not forbidden. At the point J = 2t the model (called supersymmetric) is 
associated to a graded Lie algebra and it is exactly integrable. Their integrability was 
studied in 0], although equivalent systems had been solved by other authors as in 
and 0. The ground state and the excitation spectrum were investigated in and 
low-lying excitations close to half-filling were treated in 0. The integration of the model 
by means of the nested Algebraic Bethe Ansatz (NAB A) , [|ll| in the framework of the 
Graded Quantum Inverse Scattering Method (GQISM) [|12[ was established in |]T^. The 
completeness of the Bethe states was considered in and the properties of the model in 
an external magnetic field were studied in [ p3| . 

In this paper we propose a model with two different kinds of sites. The sites of the first 
type may be unoccupied or occupied by an electron with spin up or down, but the double 
occupancy is not allowed. We will call these sites as t-J. The sites of the second type may 
be occupied by an electron with spin up or down, or by two electrons with antiparallel 
spins in a singlet state, but vacancies are not possible. We will call these sites as "J-t" 
because it is as an inverse t-J model. We can look these sites as frustrated Hubbard sites. 



Other t-J models with impurities have been proposed by different authors [|T6 | 

In order to make the inhomogeneus system we have used the i?- matrix of the t-J 
model and the Graded Yang-Baxter Equation (GYBE). We have take other solution that 
we associate to the " J-t" sites, that solution can give us another integrable homogeneus 
system. The system that we propose is formed by alternating Nh states in the t-J sites 
with Np states in the J-t sites. 

As we will see, the method is trivially dependent of the number of sites of each 
class, then we are going to diagonalized the transfer matrix for a general lattice with Nh 
and Np states in the respectives sites, although we have used an alternating chain with 
Nh = Np = N/2 to compute the Hamiltonian. This Hamiltonian is the sum of nearest 
neighbor interaction terms (two site operators) and next-to-nearest neighbor interaction 
terms (three-site operators). The diagonalization is made using the method that we pro- 
posed in ref. [0. This method is more general than the usual NABA, and it has been 
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used by other authors in a different system [18|. Recently J. Links and A. Foerster have 
proposed a model with the same kind of alternating site states that in our model and 
they solve it with the same method. 

As we said before and will see along this article, all resuls that we have, are independent 
of the number and position of the two kinds of sites, then our system can be considered as 
at- J model with impurities J-t. 



2. The model 



The t-J model may be deduced, in the framework of the graded Lie algebras, by means 
of the L-operator. We are going to take the operators given by Efiler and Korepin in ref. 

m ■ 



where 



ra,b 



^a,b5c,d 



(2.1) 

(2.2a) 
(2.26) 



and ej are the Grassmann parities of the basis vectors. We will use in this paper the FFB 
basis, that is, ei = e2 = 1 and €3 = 0. The upper indices are in the space that we call 
auxiliary, and the lower indices in the site space. 
The operator L^*"'^) verifies the GYBE 



R{X- H) 



L(*-^)(A) ®L(*-^)(;u) 



R{X-i,), 



where the i?-matrix is given by 



R{X) = XP + il, 

and the tensor product is the graded tensor product, which is defined as 



(2.3) 



(2.4) 



(2.5) 



Equation (|2.3| ) in components is. 



^(A - ^)^£^^*-"HA)-^^L(*-^)(^)f^^;^(-l)^^i(-i+-) = 



In order to make a mixed lattice, we need other L-operator associated to the J-t sites. 
That operator must fulfil the same GYBE eq. ( |2.3|) . Inspired in the same methods given 
in ref. [0 , we have found that 



L^-'-'\\) = {\--)I-iQ (2.7) 



verifies our requirements, being 



We have used the t-J basis 





|0) = , li) = 1 , I T) = , (2.9) 




and the J-t basis 





it) = , I i) = 1 , I T) = , (2.10) 




It is easy to show by direct calculation 

L(^-*)(A)^>(^-*)(-A)f;J = p{\)5^,.5^,, (2.11) 

with 

p{\) = -\~X\ (2.12) 
We build the monodromy matrix alternating the t-J and the J-t L-operators. 

T„,,(A, w) = l(*-^)(a)-x^^'-*Ha + wr^;%L^'-'\\r:;z ■ ■ ■ 
. . . l(*-^)(a)::^-^,^""^l(-^-*)(a + wY^'-'i X 

EN/2, , , v-^iV/2, , v-^iV/2-1, . , v-^iV/2-1, ,"| 

^ ,^1 (««+''0 2^^.^l(«J+«J + l) + 2^,^l 2^^-^i («i + l+ai + l)J j^2.13) 

with N even and aj = ctj = if j > N/2. We have used the standard matrix product in 
the auxiliary space and the graded tensorial product in the site space. 
In the auxiliary space we can write the monodromy matrix as 

Ai,i(A,m;) Ai,2(A,m;) Bi{X,w) 
T(A,«;)= I A2,i(A,«;) A^^^^w) B2{X,w) \ (2.14) 
Ci{X,w) C2{X,w) D{X,w) 
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The transfer matrix is given by 



F(A, w) = strace(T(A, w)) = D{X, w) - Ai,i(A, w) - A2,2{X, w), (2.15) 

and the verification of (GYBE) by the monodromy matrix assures the commutation of the 
transfer matrices for different values of the argument, that is 

[FiX,w),Fi^,w)] = 0. (2.16) 

The corresponding associate Hamiltonian is obtained by taking the first logarithmic 
derivative of the transfer matrix at equal zero spectral parameter 

HH = -iJ-^HF{\^))\x=o^ (2.17) 

where J is a constant. 

We write this Hamiltonian as the sum both operators: the nearest neighbor interaction 
term and the next-to-nearest neighbor interaction term 

. J N-l . J N-l 

^W = ^i:e. + ;^Eeu... (2.18) 

j odd j odd 

with R{0) = ciI and 

KJ+i)„ = L^'-'\^)rsL^'-'\-^t} (2.19a) 

The dot over the operator shows the derivative. 
Taking w = we find that, 

N-l . . 

J odd 

+ ^ {dj+i^^Cj+2,a - h.c.) - 2nj+inj+2 > - 6nj - 4nj+i 

+ 3{ E (4,<.^.-+2,a + h.c.) + iS-S+^, + S+S-^, + 25|5|+2)} 



(2.20) 



where and n are the creation operator and the number of electrons operator in the t-J 
sites, and and n are the same operators in the J-t sites. Also, 5' is the spin operator. 
In terms of the elements of the L*""^ operators, they are 



1 



3,1 

.b 



t 

V 

t 

-^^{L(--)(A)i;i + L(-^)(A)^;^4 + ^^^l^^ (2.21) 



X-i 

and in terms of the L"^"* 



t 
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i 






_lr(>^-*)('Al^'^ 




\dA - 






L"-^ \a,b 






Ja,& - 


ir(^-i)fAl2'l 




[sx,i, = 













2A-3i 



We have also used the following notation, 



a 



^(T)=i, r{[)=], e(T) = l, e(i) 



-1. 



(2.22) 



(2.23a) 
(2.236) 



Obviously, the Hamiltonian ( |2.2UD is not hermitian, however as we will see in section IV, 
it has real eigenvalues. In refs. and [^.is show that the supersymmetric t-J no 

hermitian hamiltonians with quantum group invariance enjoy this property. In a most 
general case, in ref. |2^] complex hamiltonians with PT invariance is proof that they have 
real spectra. 



3. Algebraic Bethe Ansatz 

The monodromy operator T verifies the GYBE, independently of the combination of 
t-J and J-t sites that we take in our system. Then, we are going to take sites of the 
first type and Np of the second. The space of states of the total system will be. 



^t-J j:pJ-t 
j 



E = ^Et' 6^EJ-^ (3.1) 



where El is the space of states of the site i and E'^ * the space of states of the site j. 
In order to diagonalize the Hamiltonian by means of the Nested Algebraic Bethe 



Ansatz (NAB A) , we should need to find a state of the system that verifies [12 



A,Jv)^5,^j\\v), (3.2) 
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but this is not the case. To overcome the problem, we must foUow a modified Nested 



Algebraic Bethe Ansatz MNABA, which is described in reference [0. For this purpose, 
in a first step, we build the vacuum subspace 

^^ = (8)|0).(8){eb, (3.3) 

i j 

where {e}j denotes the subspace generated by the vectors | [) and | t) in a J-t site j. Any 
\\w) vector in the vacuum subspace verifies 

A, ^j{X)\\w) e , ^,J = 1,2 (3.4a) 
C,{\)\\w)t^Q , z = l,2 (3.46) 

B, {\)\\w)=Q , i = l,2 (3.4c) 
D{\)\\w) = [a'{\)f' MX)]^" \\w) (3.4d) 



where G O and 



a'(A)=A + z, hi{X) = X--. (3.5) 
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From the GYBE we get the relation of commutation 

A„,,(^)a(A) = {-lr-'-g{^i - A)r(^ - A)^;^Cp(A)A„,,(^) 

+ h{^-\)CMAa,c{\). (3.6a) 

D{^i)C,{\) = g{\ - iJi)C,{\)D{iJi) - h{\ - ii)C,{iJi)D{\), (3.66) 

C^{\)CM = r(A - ^) J^a(^)C,(A), (3.6c) 



where 



g{^) = (3.7a) 
hill) = -, (3.76) 



In order to solve the equation 



F(A)^' = A(A)^', 
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(3.8) 



we build the state, 



vl/(A) = vl/(Ai, . . . , A,) = CaA^i) ■ ■ ■ C„^(A,)X"i-"'^||l), 



(3.9) 



with II 1) G O. When we apply the -D(/u) operator to using ( pl)a — c\ ), we push this 
operator to the right of the C operators, and we get a wanted term characterized by the 
C operators conserve their arguments, and several unwanted terms characterized by the 
arguments of the C operators are interchanged. The wanted term is. 



j=i 



and the k-th unwanted term is, 



(3.10) 



...a,_,(A,_i)M('=-i)(Afc_i)^-,tX"^'---'"1|l), (3.11) 

with M given in the appendix A. 

The application of „ to the state \E' is a little larger but straightforward. We get 
again wanted and unwanted terms, and after some calculations we find the wanted term. 



n gil^ - Xj)C,, (Ai) ...Cp^ iXr)A,M [Z(^, A)^-;.::;^: 



X 



b.a 



1). (3.12) 



where the Z operator is 



= L(i)(^ - XrY.^-^L'^'Hfi - Xr-i)iz\:t--l ■ ■ 



(3.13) 



with 



and P(i)};;; = 5e,Aj(-l)^'^'^=. 

Due to J = 1 or J = 2, we have = 1, then ( p.l4|) becomes 



(3.14) 



^^'^(A)ff = -r(A){, 



b 

e,d 



(3.15) 
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Now we will get ready for prepare the second step of the MNABA. We define the 
monodromy matrix at second level as 

T(2)(^,A) = A(^)-Z(^,A). (3.16) 

In the auxiliary space, that now his two-dimensional, this operator can be written as 

The transfer matrix at second level is 

f(2)(^^X) = strace [T(2)(;U,A)j = -A^^\fx,X) - D^^\^,X). (3.18) 
Now, the wanted term ( p.l2|) can be written as, 

r 

and the k-th unwanted term as 

r 

Ki^ - Afc) n 3^^^ ~ ^j)CpMCpu+A^k+i) ■ . .Cp^{\r)Cp^{\i) . . . 

i = i 

...Cp,_,{\j..,)M^'-'\\,,.,)ll-^^^^ (3.20) 

We have found a new problem; we have to solve the equation, 

F(2)(^,A)X||1) = A(2)(^,A)X||1). (3.21) 

For this purpose we use the fact that L*^^)(A) verifies the GYBE with the r-matrix, and 
then we can solve ( p.21|) following the same path applied before. 
The new vacuum is 

l|l') = (8)|0i).(8)(j), (3.22) 

where the state |0 J. ) is composed by a state |0) in the t-J sites and | | ) in the J-t sites. 
The first tensorial product is the space where A{^) works, whereas Z(//, A) works in the 
second one. This vacuum verifies 

S(2)(^^A)||lO = , (3.23a) 
C(2)(^,A)||1') ^ , (3.236) 

A(2)(^,A)||lO = [%)]^M^i(^r^n4^S^I|l') , (3.23c) 

f^^ ^(^ - Aj) 



with 



ai(A) = A+- 



6(A) = A. 

The new relations of commutation are obtained from the GYBE for T^^) 



D(2)(^)C(2)(A) = g{X - ^)C^^\\)D^^\^) + - \)C^''\^)D^''\\), 
In the second level, we build the state 

Applying A'^^^ to this state we get a wanted term 

i=i j=i 

and unwanted terms of which, the /c-th is 



xC(2)(^)C(2)(^,+i)...C(2)(^,_i)||l') 

Also for 

we get a wanted and unwanted, the wanted is, 

i=l ' 

and the /c-th unwanted term is, 



f (/^ - Aj 



.Vl>(2). 



1 



9{^^k - Aj) 



i=l 



(3.24) 

(3.25a) 
(3.256) 

(3.26) 
(3.27) 



(3.28) 



(3.29) 



(3.30) 



The cancellation of the unwanted terms, (|Tl|) with and ( ^^281) with ( |O0|) , 

gives us the ansatz equations 



6(Afc) 
a'{\k) 



No 



1 Nr, 



n 



\ fi'(Afc - in 



/C = 1, 



n 



gilin - in) 



n = 1, . . . , s. 



(3.31a) 
(3.316) 
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On the other hand, collecting the wanted terms we get the eigenvalue of the transfer matrix, 

r 



[bMf^ n^(^ - n 9i^i - ^) + Il9(^'^ - ^) r (3.32) 

i=l j=l i=l 



4. Thermodynamics of the model 

The eigenstates of the transfer matrix can be characterized by observables operators 
that commute with it, and then with the Hamiltonian. We have found that the following 
observables 

7V-1 

Oi = J]](ni + ni+i), (4.1a) 

i=l 
i odd 



i=l 
i odd 

N-l 



i=l 
i odd 

N-l 



i=l 
i odd 



commute with the transfer matrix 

[O,,F(A)] = 0, z = l,...,4. (4.2) 

We must note that Oi is the total number of the electrons, O2 is the difference between 
the number of holes and number of electron pairs, O3 is the number of electrons with the 
spin up in the {|0 >,| |>,| t>} base plus the number of electrons with the spin down 
in the {| i,t>,| i>,| T>} base and O4 is the complimentary of O3. Only two of these 
quantities are independent, the others are linearly dependents and thus, we are going to 
use Oi and O3. 
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The commutation relations of the Ci creation operators with these observables aUows 
us to relation their eigenvalues with the numbers r and s of the ansatz equations. They 
are: 



[Oi,C,] 


= 






[03,C,] 




[04,C,] 





(4.3a) 
(4.36) 
(4.3c) 
(4.3rf) 



and a long but straighforward calculation, gives us the action of these operators on the 
eigenstate ij{Xi, ■ ■ ■ , A^), 



Oi^(Ai,- 


■ ■ 5 A^) 


= {r + Np)i;{Xi,---,Xr) 


02^(Ai,- 


■■,A.) 


= {~r + NhmXi,---,Xr) 


03^(Ai,- 


■ ■ 5 A^) 


= {r-s + Np)ij{Xi,---,Xr) 


04V^(Ai,- 


■■,A.) 


= sV^(Ai, ■ ■ ■ , A,.) 



(4.4a) 
(4.46) 
(4.4c) 
(4.4d) 



The next step is to solve the equations of the ansatz in the thermodynamic limit, 
i.e., when A^, A^^ and Np are infinity, but Nh/N and Np/N remain finite. For this, it is 
convenient to reparametrize the roots as follows 



(2) 



I 

2' 



I,---, 



(4.5a) 
(4.56) 



and then the equations (|3.31a, 6| ) are written 



r4" 




i " 
2 


Nh 


k> 


+ 


i 

2 ^ 






+ 


i " 
2 








i 

2 ^ 





n 



,(1) 



'^J 2 



,,(1) _ ,,(2) , 1 



n 



,(2) 



,(2) 



n 



,(2) 



(1) + 1 J-l ..(2) _ (2) _ . 



The energy is obtained with 



d 



E = -zJ— InA(A) 
uA 



-iJ 



A=0 



r 



a'(0) 



^>i(0) ^tt 



(4.6a) 
(4.66) 



(4.7) 
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and using (4.5 ) 



It is convenient to introduce the function 



[X) 



2 arctan (x) = —i In 



1 + ix 
1 — ix 



and taking logarithms in the equations of the ansatz ( [4.6a, we can write 



r(l) 



(2) 



,(2) 



27rl 



(2) 



(4.8) 

(4.9) 

(4.10o) 
(4.106) 



with and I^i^'* integers or half-odd integers and each set of these numbers determines a 



solution for roots v^}^ and Vf^ . 

In solving a BAE set, we have real roots and complex conjugate roots grouped in 
clusters that we call n-strings. The solution for the ground state in the thermodynamic 
limit can be obtained minimizing the free energy distribution . In the Heisenberg model 
with arbitrary spin is found that the solution for the ground state is formed by two-string 
roots [^. In the t-J model, the two level roots are mixed, they proliferate rapidly and it 
become difficult to determinate the roots which parametrize the ground state. Numerical 
analysis suggest that the structure of solutions in the ground state are two-string roots . 
We assume the same hypothesis in our model. 

Then, the solution in ground state is given by two-string roots and for the lower excited 
states some roots are real. The equations ( [4.1Ua, fc| ) can be reparametrized supposing that. 



(2) 



in the r roots of the first level, v\^\ j = 1 



V 



(1) 



r, there are 2ri roots of the form 
/ = l,---,ri (4.11) 



with €i real, that are grouped in ri two-strings. The rest r — 2ri roots are taken real and 



they are designed by. 



V 



(1) 



1 = 1, 

(2) 



(r — 2ri). 



(4.12) 



With the same assumptions, the s roots v) of second level are written 



(2) 



A ±2' / = l,---,si 

ui. l = l,---,{s-2si] 
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(4.13) 



with (3i and vi both real. 

Taking this in account, the equations (|4.10a-5|) result for the respective roots. 



s — 2si 



Sl 



(1) 



(4.14a) 



1=1 

s — 2si 



1=1 



Sl 



NH(piek) - J2 ^(^^ -^^)-Yl ["^(^efc - 2A) + (P{-ek - -A)] = 27r jj^'^ (4.146) 

i=l 1=1 

r— 2ri ri s — 2si 

Np(l){2vk) + ^ (j){2i'k - 2rij) + Y(j){i'k - ej) - ^ (j){uk - Ui)- 



Sl 



i=l 



[(/>(2z/fc - 2A) + </>(T7^^fc - ttA)] = 27rli- 



(2) 



(4.14c) 



z=l 
r — 2ri 



ri 



s-2s 



Y m2/3k-2ly,)+<Pi^Pk-lly^)]- 

i=l 

Y mf^k - A) + H^iPk - A))] = 2nj 

i=l 



(2) 

k 



(4.14d) 



where we have used the identities in appendix B. 

In the thermodynamic limit the system is described in the language of distribution 
functions of roots in both levels with particles and holes. In our case we define it as follows: 



pi = lim — ^ 

N^oo N{r]j+i-r]j^ 

P2 = lim 



0"i = lim 



1 



02 = lim 



iv^oo A^(ej_|_i — ej) 
1 



(4.15) 



and we call [-Bi], [-62], [C'l] and [C2] the regions where pi and p2, (Ti and a2 are defined 
respectively. 

As usual[|25|, we introduce the functions, 



ZpiiX) 
Z.iiX) 



1 



2txN 
1 

2'kN 



s — 2si 



i=l 
s — 2si 



1=1 



Sl 



NhHX) - 5^ (/>(A - u,) - J2 [H'^X - 2A) + <P{-X - -(3, 



i=l 



i=l 



(4.16a) 
(4.166) 
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27riV 



r — 2ri 



s — 2si 



j=i i=i 



Sl 



5^[(/>(2A-2A) + <^(|a-|a 



27rA^ 



3 3 

r—2r\ r\ 



(4.16c) 



+ ^ ^(A - ry,) + [./>(2A - 2e,) + 0(-A - -e,)]- 



[</>(2A - 2Z.0 + </>(|a - |z.,)] - £ [2<^(A - A) + (/>(^(A - A))] 



1=1 



(4.16d) 



These functions are monotonyally increasing in their respective definition regions 
[S1/2] and [C1/2], and their values are integers when A takes the value of a root. Be- 
sides, there are other values of A where the Z functions take a integer or a half-odd integer 
value, but they do not correspond to a root. We call each of these values a hole. The 



distribution functions are the derivatives of the Z functions 25 



(4.17) 



where 6i are the position of the corresponding holes. The thermodynamic limit is obtained 
by doing, 

^ / ^A/(A)p(A). (4.18) 



Then, from (|4.16a-d|) , the corresponding distribution functions in this limit are. 



N 



Pi(A) 



1 

2^ 



2N f 

2cP'{2X)-— / cP'{2X-2^)p2{fi)dfi 



1 

2^ 



— / (p'{X- fi)a2{n)d^ 

^^h J[C2] J 

N f , 

2</>'(A) - — / (f)'{X- fi)p2{p)dp 

J-^h J[B2] 

2N r 12 
- — / [<^'(2(A - p)) + -<^'(-(A - p)]a2{p) 

J^h J[C2\ 

N If 2N f 

— P2(A) = — 2</>'(2A) + — / </>'(2A - 2p)pMdli 
N f N f 



(4.19a) 



(4.1%) 
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2N 



P J[C2] 



N 



a2(A) 



1 

2^ 



N 



/,'(2(A-^)) + i0'(^(A-^)]a2(^) d^i 



(4.19c) 



+ 



2N 



P J[Ci] 



P J[Bi] 



2N r 1 2 

+ — / [<^'(2(A - ^)) + -c^'i-iX - ^)]p2(^) dp 



[2<^'(A-^) + ^</>'(^(A-^))]a2(^) dp 



(4.19d) 



Quantities corresponding to observables defined before, can be found on this state in 
the thermodynamic Umit in the same form, so the energy is given by 



E 
N 



J 



-Nh + 2iVp 
N 



+ 2 / (f)'{2ij.)pi{p) dp + / (f)' {p)ai{p) dp 



The Oi index (number of electrons) defined in (|4.1a| ) is given by 



He. =-- 



Oi 

N 



r + Np 
N 



^ + I pi{p) dp + 



N 



(Ji{p) dp. 



(4.20) 



(4.21) 



The difference of magnetization 5"^ between the t-J sites minus the magnetization in the 
J-t sites (|4.1c, cj ). 



QZ QZ 

N 



O3-O4 
2N 



(r + 



2s) 



= ^-/ p2ip)dp-2 a2{p) dp.. 
^ J{Bi\ J{Ci\ 



Using ( |4.196, c| ) and ( [1.2U| ), we find an important relation. 



N 



Nh - 2Np 
N 



27r((7i(0)-p2(0)) 



(4.22) 



(4.23) 



In our model, as we said before, we are going to suppose that the configuration of the 
ground state is given by two-string roots in both levels of the BAE, as it is suggested for 
the t-J model in ref. ||]. With this hypothesis, since pi = p2 = 0, the equations ( |4.19a-d| ) 
are reduced to two equations and are simplified considerably. 



JNh 2tc 



(4.24a) 



Tr^2(A) = — 
Nh 2tc 



2N r 12 
2</''(A) - — / W{2{X - p)) + -c/>'{-{X - p)]a2{p) dp 

J-^h J[C2\ ^ ^ 

<P'W + / [<^'(2(A - ^)) + tt</''(tt(A - p)]aM dp 

NT 11 

- — / [2cP'iX -p) + -<P'i-{X - p))]a2ip) dp (4.246) 
iVp J[Ci] ^ ^ 
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and the integration regions can be determined by Ue and s^. 



N 

ne = -rf + 2/ aMd^. (4.25a) 

■'[C'l] 



^-2/ a2{fj.)di,. (4.256) 
^ ■'[C2] 



These equations can be solved numerically for any filling rig and magnetization s^. 
If we take a system where the integration limits turn out to be, 

[C^i] = [C^2] = (-00,00), 

our equations can be solved analytically using the Fourier transform and some identities 
in Appendix B. 



Let be 



then, we obtain. 



/oo 
(Tj(a)e*"^c/a. (4.26) 
-00 



^2 (a) = ^ ^ — + ^ 3 (4.276) 

N 4cosh2 f A 8cosh^ f 

These expressions can be used to determinate the main parameters of our system 

^ = 2 J a,{a)da = 2a,{0) = ^^^^^ (4.28a) 
^ = 2 I a2ia)da = 2a,{0) = ^^^^^ (4.286) 



and using ( |4.25a-6D we find 



- (4.29a) 

= 0, (4.296) 

We must note that for Np = 0, we have rie = 1 and = 0, this constitutes a antiferro- 
magnetic state with half-filling. 

From ( ^.20| ), the energy is given by 

- = jr^(l - 21n2) + ^(^ - 3)1 (4.30) 
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and, for A''^ = we obtain the result of the t-J model. 
As we know, it must verify 



— - < rip < 

N - N 



(4.31) 



then, from (|4.29a|) , this solution will only be true if. 



Nh - 



(4.32) 



5. The excitation spectrum 



We start again from the BAE ( [4.6a, d^ . The more general string solutions will be of 
the form, 

= ^fci + m = -M,---,M (5.1a) 



^Um) = + ^ = -M'^ ■■■^M' (5.16) 



with M and M' integer or half-odd integer. 

Following the same method as before, we multiply the equations of the same string, 

(i) 

and we obtain that their center m'' ^ ~ 2' verifies the equations 

2Nh arctan — ^ = 2nll]l + ^^A^^" ^^^Im - ^j'm") (5-2a) 

2 M" j=l 

(2) 

,,(2) ^L' 



- 2Np arctan ' = -2nP^lj " ^j-m') 

"^2 M' j=l 



M" 1=1 



where is the number of strings in the i level, and the $ and \E' functions are defined in 



appendix B. A solution of ( p. 2a, 6|) , is determinate by specifying sets of integers or half-odd 

r(l) 1 
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integers {ij^M} and the regions where the {C} roots are distributed. 



As we did before, we define the functions, 

,(2) 



^mW = — arctan , ~ ^ 5Z 5Z ^A'/,m"(A - v^'^]^,,) (5.3a) 



^ 2 M" j = l 



(2) 
M' 



N A 1 

(A) = arctan — " ^ E E " 



,(2) ^ 
j,M', 



M 



-I IVl" 

+ ^EE*^M-^"(^-4m"), (5.36) 

M" j = l 

that are monotonic increasing and reach integer or half-odd integers values when A takes 
the value of a root. 

By counting the number of roots, and calling H^^j the number or holes in the sea of 
M-strings at level i, we have, 

and supposing that the centers of the strings are distributed along the real numbers, 

2^L,m + 1 = ^mM (5-5) 
and then, using ( |5.3a, (|5.4|) and (|5.5|) , we obtain, 

+H^M^ =Nu-2 J2 KiM, M")y^^], (5.6a) 

M">0 



where 



and 



j(M Ml-/2^i + I ifMi = M2, , . 

•^^^^'^^^-\2min(Mi,M2) + l if ^ , ^^'^^ 

KIM M\ /^2 + | ifM2 + i<Mi, 

Besides, the number of roots, obviously, must verify , 

5^(2M + l)i.W (5.9a) 

M>0 

s= (2M + l)zy£) (5.9a) 

M>0 



r = 

M>0 
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and now, we can apply (|5.6a, t\ ) and (|5.9a, tjj ) to the ground state, according to which type 



of string we considered that forms it. 

If we suppose, as before, that the ground state is formed with two sea of two-strings, 
we have, 



= AT, - - .J^) - 2 V .gl (5.10a) 

2 2 



M">i 



iff = TV, - - 4 Y: -'mI - -i^'' - 2 E -S" (5-10^) 



M'>i M">i 



and then, with this results, we obtain from (|4.28a, l\ ) 



ffW ^ ^(2) _ ^g -^^^ 

2 2 

that is to say, we have two two-string seas without holes. 

Under these hypothesis we can analyze different types of excitations by introducing 
holes in the ground state and keeping constant some observables. The one that we are 
going to consider is maintaining constant the electron number and the magnetization 

He = const., Sz = const. , (5-12) 

then, by imposing (|5.9a, fc|) and ( |5.10a, 6| ), we obtain that the state must be characterized 
having two real roots and a hole in level (1) 

z.^') = 2, h[''^ = 1, (5.13) 

2 

and besides, 

^(2) _ Q (1) _ 2N^-N, _ ^ (1) _ Q 

(5.14) 

o-(2) _ n ,,(2) _ 2N^+N^ (2) _ ^ 

2 2 ^ 2 

These conditions correspond to one the two states, the first one is obtain from the ground 
state by changing two consecutive sites with spin up to a new state with the two sites in 
spin down, the second one is a state with one of the two sides without electrons and the 
other with a pair. 
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Under our hypothesis about the ground state, we can calculate the contribution of 
every hole and real root in both levels, to the energy of a excited state compared with the 
energy of the ground state. Using ([4.17| ), the equations (|4.24a, must change to. 



IN C 12 

J^h J[C2\ ^ ^ 

^^2(A) + ^ E ^(^ - ^^^) = i k(^) + ]^ E '^'(^ - 

/i2 P ri 

2N f 1 2 

+ — / [</>'(2(A - y)) + -</>'(-(A - ij)]aM 

^> -'[Ci] <3 



(5.15a) 



2 



N 

n: 



(5.156) 



where 9hi parametrizes the holes at level i and 9^ the real roots. 

The system can be solved as before, by using the Fourier transform, and we obtain 
that the distribution functions can be written as. 



o-i(a) = ai(o)(a) + ai(„)(a) 
0-2(0) = 0-2(0) (a) + <3-2(n)(a) 



with 



5"i(o)(a) 
5"2(o)(a) 

^l(n)(«) 



Np e 



-3|a|/2 



1 + e-l«l (1 + e-l"l)2' 

'h e" 



g-3|a|/2 g-|a 



+ 



(1 + e-l"l)2 (1 + e-l«l)3 



(5.16a) 
(5.166) 

(5.17a) 
(5.176) 



1 

iV 



E 



.1/2 



(l + e-l°l)2' 



-E 



(l + e-i«l) 



-3|a|/2 



7-1 



(l + e-l«l)2' 



(5.17c) 



<5"2(n)(a) 



1 



E 
+ E 



|a|/2 



(l + e-l«l)2 



(l + e-l°l)3 



-E 
+E 



(l + e-l"l)3 



.1/2 



^•2 



(l + e-l^l) 



(5.17d) 
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The contribution to the energy per site is, 

From this expression, we observe that the holes in both levels and only the real roots of the 
first level give contribution to the energy. A straighforward calculation show the following 
rules: 

i. - every hole in the first level gives a contribution Ae = —ei{9hi) 

ii. - every hole in the second level gives a contribution Ae = £2(^/12) 

iii. - every real root in the first level gives a contribution 

Ae = sech(7r6'^J + e2{0r^) 

iv. - every real root in the second level gives a contribution Ae = 0, 
being 

da. (5.19) 



nN Jo et cosh'' f 

We can apply these rules to the first example, where we have held constant the number 
of electrons and magnetization respect to the ground state and it is characterized by two 
real roots and a hole. We can suppose that the roots and the hole are parametrized by the 
same value 6*, it is said, the state is coming by changing a two-string of the ground state 
by two real roots, which leave a hole in the two-string sea. Doing the integrals in Ae, we 
obtain for the state, 



-4 



^ +2Re[*(J -fz^)] -2Re[^'(^ -f4)] 



cosh(7r0) ^ '2 ' M 2- 

-4^(lm[^'(i-F^^)]-Im[v&(l + z^)]) 

-2Re[^'(A -f2Re[^'(z^)] -2(l + 21n2) 



(5.20) 



where \E'(a:) is the derivative of the logarithm of the Euler gamma function. 

Taking an alternating chain = Np = N/2, we can represent the energy given by 
( p.20| ). The results are shown in fig.l. As we can see, Ae is null out of an interval around 
the origin. Then a conclusion is that there is not energy gap between the excited and the 
ground state. 
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Appendix A. 

Using tlie relation of commutation ( |3.6c| ) we have 



C'ai(Al)Ca2(A2) . . . Ca^{Xr) = Cfoj (A2)Cb3 (A3) . . . 

. . . a.(A.)a, (Ai)G(Ai, X)^::± , (A.l) 

with 

G{u, X)^::± = r{u - Xi)l\^::r{u - A^)-'- ...r{u- KfCl-. ■ (^-2) 

Taking 

m(^)(A,) = ^(A,, A) ■ ^(A,-!, A) ■ . . . ■ ^(Ai, A) , (A.3) 

we have the relation 

C'ai(Al)Ca2(A2) . . .Ca^{\r) = Cb^{Xk) ■ ■ ■ Cb^{Xr)Cb-^{Xl) ■ ■ ■ 

. . . Cb,_, (A.-OM^'^-i) (Afc_i)^^';;;i; . (A.4) 

Appendix B. 

In this appendix, we are going to give brief remarks about the function, 

(pix) = 2arctan(a;) = — zin , (B.l) 

1 — ix 

that takes values in the interval —n to +7r when x goes from — oo to oo. 

The function has the following properties, that can be proof by a straighforward 
calculation: 

(j){x + i) + (j){x - i) = TV + (B.2) 

(/)(x + 2z) +(^(x-2i) = (^(-) -(^(x) (B.3) 

3 

</>(x + ^) + </>(x - ^) = '/'(f ) + ^i"^^) (B-4) 
The Fourier transform that we have used 



1 



oo 



p(A) = — / p{a)e'''^da, p{a) = / p{X)e-'"^dX, (B.5) 
27r ./_oo 
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The derivative of function (j) is 

d(t){X) 2 



dx 1 + A2 

and then, the Fourier transform, 



It is convenient to define the functions 

M2+M1 



*Mi,M2(a;) = 2 ^ [arctan — h arctan — -j-y] (B.8) 



X X 

[ arctan 

n=|M2-Mi| 

and 

Ml . 



n=— Ml 



(B.6) 



/oo 
0'(A)e-^«^dA = e-l"l (B.7) 
-00 



X T III' 

$Mi,M2(a;) = 2 2^ arctan— —y (B.9) 

M2 + 2 
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Fig. 1 

Energy of the Excited State 
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